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The statistical mechanics of directed line-like objects, such as directed polymers in 



an external field, strands of dipoles in both ferro- and electrorheological fluids, and flux 
lines in high-T c superconductors bears a close resemblance to the quantum mechanics of 
bosons in 2 + 1 dimensions. We show that single component and binary mixture critical 
phenomena in these systems are in the universality class of three dimensional uniaxial 
dipolar ferromagnets and ferroelectrics. Our results also apply to films of two superfluid 
species undergoing phase separation well below their A-points near T = 0. In the case of 
directed polymers and electrorheological fluids we analyze the effects of free ends occurring 
in the sample as well as a novel directionally-dependent compressibility 
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1. Introduction and Summary 

There has recently been renewed interest in the analogy between the classical sta- 
tistical mechanics of lines and the quantum mechanics of bosons in 2 + 1 dimensions. [|I| 
Configurations of directed lines can be mapped onto the world lines of a collection of bosons 
in one fewer dimension. Examples include flux lines in high temperature superconductors 
P], polymer nematics in a strong external field, and ferro- and electrorheological fluids 



These systems can be modeled by a free energy of the form 
dzd 2 r 
The partition function 
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is subject to the constraint 



Z= f[dt][d5p]e- F/kBT (1.2) 



d,5p + Vrt = 0, (1.3) 



where Sp(r, z) = p(r, z) — po is a cross-sectional areal density deviation from its average 
value po in a plane perpendicular to z, and t(r, z) is a tangent field in the same plane. 
Given a microscopic description of directed polymers, via their trajectories defined by 
single-valued functions {rj(z)} along the z-axis, p(r,z) and t(r,z) are defined by (see 
Figure 1) 

p(r,z)=^r5 2 [r- ri (z)] (1.4) 

i 

and 

t(r,,) = £^^ 2 [r- ri (z)]. (1.5) 

i 

We assume that, as external parameters such as field strength, osmotic pressure, 
temperature, solvent quality, etc. are varied, the polynomial part of (|1.1|) in dp exhibits 
a line of first order phase transitions from a high to low density phase, terminating in a 
second order phase transition. Along the critical isochore, w = 0, and b then changes sign 
at the mean field critical point. Here and throughout, we distinguish the z direction as the 
direction along which the lines are aligned, and denote by V x the derivative perpendicular 
to the z-axis. Additionally, we use boldface to denote vectors in the perpendicular plane. 
In the case of quantum mechanical bosons, we think of z as a time-like axis along which 
particle world lines are extended. Planck's constant % plays the role of "temperature" in 
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this system, f3% is the sample thickness, and the coupling h is the inverse "mass density" 
of the bosons. 

The constraint (|1.3| ) in the context of polymers corresponds to the absence of polymer 
free ends [||],|[] , that is, conservation of polymer number as we would move up the z-axis. 
In the context of quantum mechanical bosons, this constraint merely sets t to be the 
canonical momentum density of the particles ||. Models similar to (LI) were first posed 
as Landau-Ginzburg type theories 0||j8[] but were later shown to be a direct consequence 
of the boson mapping || . 

In fact, this model is equivalent to a three dimensional uniaxial dipolar ferroelectric 
(or ferromagnet) near its critical point, first studied in 1969 by Larkin and Khmel'nitskii 
M. To see this, we use the constraint ( p..3| ) to solve for the longitudinal part of t. Upon 
integrating out the transverse part of t we obtain an effective free energy 

F ^=lJ^Wf [hi + e q l + a q l + bySp( q± , qz )\ 2 + lJdzd^(Sp)\ (1.6) 
where the quadratic part has been written in Fourier space. Here and henceforth, we 



restrict our attention to the critical isochore, and set w — 0. Equation ( |1.6| ) is the familiar 
form of the Landau theory for a uniaxial Ising ferromagnet or ferroelectric with a mean 
field critical point at b = || . 

A model very similar to ( |1 . 1| ) applies to polymer nematics in which the external field 
is turned off, so that the alignment represents spontaneously broken symmetry instead 
of being externally imposed H@|,|||| . In this case, t(r,z) represents a nematic director 
field and the term \ht 2 is replaced by the usual nematic gradient free energy []IIJ. Simple 



power counting arguments for this model suggest that nonclassical critical behavior results 
below (3 + l)-dimensions. A complete treatment, however, would require a discussion of 
the effects of hairpin turns in the polymer configurations 0. We restrict our attention here 
to polymer nematics in a strong external field, and to length scales less than the spacing 
between hairpins, to avoid these complications. 

1.1. Electrorheological Fluids and the Effect of Free Ends 

The effective attraction for polymers in a nematic solvent is similar to the fluctuation 
induced attraction of strings of colloidal spheres discussed by Halsey and Toor |TT|] . This 



attractive part of the interchain potential in electrorheological fluids can be offset by a 
hard-core repulsion at short distances, and may be adjusted by varying the dielectric 



constant of the solvent. Since an electric field is required to create directed strings of 
dielectric spheres, hairpin configurations will be greatly suppressed in electrorheological 
fluids, as there is no longer an up-down symmetry. We then expect that ( |1 . 1| ) is especially 
apt for these systems, once free ends are taken into account. Halsey and Toor [|12|, by 
considering the bulk energetics of the system, have argued that at sufficiently large fields, 
the electrorheological strings will phase separate, leading to two-phase coexistence. The 
critical point discussed here would then describe a transition to a uniform density. 

To study this point further, consider the free energy at fixed potential due to the 
electric field 

U = -— [ dVe[n(x)]E 2 (x) (1.7) 

8% J 

where n is the full three-dimensional colloid density, and e(n) is the effective dielectric 
constant. There is an overall minus sign because this is the energy at fixed E. Since 



e"(n) > |13|], it is energetically favorable to phase separate into dense and dilute regions. 
However, at finite temperature, configurational entropy plays a role. Certainly, at fixed 
electric field, the temperature can be increased until entropy dominates and there is only 
one uniform phase of directed lines. However, entropy also favors shorter chains at such 
temperatures, which could obscure the entanglement effects which are the main focus of 
this paper (see below). 



Recent numerical simulations by Tao [14] suggest that at fixed temperature there are 



actually two critical electric fields, E C1 and E C2 (E C2 < E C1 ). Below E C2 the fluid consists 
primarily of a gas of monomers. Above E C1 the colloidal particles clump together and 
phase separate into "columns" and dilute chains as Halsey and Toor's argument would 
suggest. The columns may in fact be crystalline []T4]JT5|1 . Between E C2 and E C1 there are 



dense strings of colloidal spheres with random positions in any constant-^ cross section. 
The simulation suggests a first order phase transition from the random line phase to the 
clumped phase. At lower densities, experimental measurements do not suggest crystalline 



order in the columns [16]. Thus, at low enough densities we can have coexistence of 
liquid-like columns and a gas of chains, possibly terminating in a critical point. 

One shortcoming of our model in its present form is that the chains it describes must 
start and stop at the boundaries of the sample. Both directed polymers and colloidal 
chains should be allowed to have free ends within the sample. As discussed in Section 4, 



free ends may be incorporated by relaxing the constraint ([O]) and using the free energy 



F = / dzd 2 x 

\tf + e - (v,m 2 + f (<^p) 2 + ^ P 2 + ^/ + |(a z 5 P + v x -t) 2 

By taking G — > oo, we recover our original model. Additionally, it can be shown that 
G ~ t~ x as I — > oo, where £ is the chain length ||,|J. The crossover from elongated objects 
to point objects when the chains are finite will be discussed in detail in Section 4. 

1.2. Phase Separation in the Flux Liquid Phase 

Flux lines in superconductors may also be described by an analogy with bosons ||. 
However, in the absence of magnetic monopoles, flux lines cannot have free ends inside the 
sample. Their behavior, then, is much closer to the behavior quantum mechanical bosons, 
as the constraint (|1.3j ) is always satisfied. 

Calculations to date for fluctuating flux lines have assumed a purely repulsive pair 
potential [0 . Muzikar and Pethicke WJ\ , however, have argued that an effective attractive 



interaction arises under some circumstances near H C1 in the extreme type II limit. This 
interaction would convert the continuous onset of flux penetration at H C1 to a first order 
transition. An even more intriguing situation would arise if the interaction between the 
flux lines had two distinct minima (see the inset to Figure 3). In this case, there should 
be a first order transition at H C1 , and a line separating two coexisting flux liquid phases, 
as in Figure 2. Our theory applies to the critical point marking the terminus of this first 
order phase boundary. The critical isochore is shown as a dotted line in the figure. 

1.3. Quantum Critical Phenomena 

Our analysis of course, can also be applied directly to boson systems at T = 0, which 
are equivalent to infinitely long directed polymers |I| ■ In the case of a single particle species, 
we would be describing the terminus of a line of first order transitions separating superfluids 
coexisting at two different densities. How this situation might come about is illustrated in 
Figure 3. The inset shows an unusual pair potential with two distinct minima describing 
particles which, like helium, are light enough to avoid crystallization at T = at low 
pressures. The two minima lead to two distinct liquid phases with increasing pressure. At 



low pressures, the stable zero temperature liquid has near neighbors occupying the deeper 
outermost minimum. As pressure is increased, there is a first order phase transition to a 
denser liquid with the innermost minimum occupied, followed eventually by crystallization 
at even higher pressure. The two liquids must be superfluids at T = 0. In Figure 3, the 
parameters have been adjusted so that the critical point for the dilute and dense superfluids 
SI and S2 occurs below the dashed A-line describing the transition to a normal liquid with 
increasing temperature. There is a second critical point associated with the coexistence of 
this normal liquid with a dilute gas. 

To obtain the critical point of interest to us here, we must increase % (or decrease the 
mass) at T = until the two superfluid liquid phases become indistinguishable as at the 
point c in the figure. Although this gedanken experiment is clearly rather esoteric, some 
control over quantum fluctuations may be achievable for two-dimensional helium films, 
where out of plane fluctuations alter the effective two-body potential. By varying the type 
of substrate, it may be possible to create a purely repulsive effective potential and a novel 
superfluid gas phase at T = [ITS]. We should stress, however, that we do not expect the 



unusual phase diagram of Figure 3 for either purely repulsive potentials or conventional 
helium potentials with a single minimum. 

The density operator of the bosons described above is given by ( |1.4|) . Upon differen- 
tiating with respect to z we have 

a z p(r, z) = - Y, ^P- ■ V.,r5 2 [r - r,(*)] = -V x • £ ^5 2 [r - r,(z)] (1.9) 



dz ' " t—^ dz 



i 



The final sum on the right hand side is just (|1.5| ), so we have derived ( |1.3| ). The quantity 
t(r, z) = Y^i dz " [ r — r i( z )] IS J us t the momentum operator of the bosons [HH]. On a 
more formal level, note that the "time dependence" (i.e. the z dependence) of the quantum 
operator ( |L4] ) is given by the (Wick rotated) Heisenberg equation of motion, 

d z p(r,z) = [p(r,z),H], (1.10) 



where Tt is the underlying quantum Hamiltonian. Upon identifying ( |1 . 1| ) with the integral 
over imaginary time of the coarse-grained hydro dynamic Lagrangian associated with this 
Hamiltonian (as is appropriate for finite temperature quantum statistical mechanics), we 
recover ( p..3| ) from Ql.lU ) upon using the commutation relation between p(r, z) and t(r, z) 



[p(r, z), t(r', z)] = ^V,5 2 (r - r'). (1.11) 
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In the hydrodynamic limit considered here, this commutator becomes the classical "Poisson 
bracket" relation between density and momentum, which arises because the momentum 
operator is the generator of translations | I5[| . The constraint ( |1.3| ) thus directly reflects 



the non-commutivity of the underlying position and momentum operators. 

Note that our results represent a new universality class for Ising-like quantum critical 
phenomena. A related type of quantum critical phenomena, extensively studied in the past 
|20[1 , arises for the Ising model in a transverse magnetic field. The quantum Hamiltonian, 



defined on a hypercubical lattice in (i-dimensions is 






n = -J^ctc]-h^c* (1.12) 

(M> 3 

where the {a^} are Pauli spin matrices. The ordered state at T = is destroyed for 
sufficiently large h. This disordering transition is known to be in the universality class of 
the (d + 1) -dimensional classical Ising model POP - Equation (|1.12| ) describes a quantum 
lattice gas, with the world lines of the "particles" (up spins, say) in a T = path integral 
formulation representing the directed polymers discussed in this paper. Because the cxj is 
a spin flip operator, these "polymers" stop and start at random in (<i+l)-dimensions. The 
appropriate coarse-grained Landau description would be similar to (|1.1|) without, however, 
the constraint ( |1.3|) . The spin flip operator in (|1.12| ) does not respect this conservation law, 
so it is not surprising that we find a different universality class for boson critical phenomena. 
Introducing free ends in our model has the same effect as the spin flip operator in ( |1.12|) . 
As discussed above, we then find a crossover from the (d+ 1) -dimensional uniaxial dipolar 
Ising critical phenomena to the (d + 1) -dimensional Ising behavior expected for ( |1.12|) . 

1.4- Directed Polymer Blends and Quantum Binary Mixtures 

By extending the original model to allow for binary mixtures, we can study the con- 
solute point of two different directed polymer species. In particular, we can determine 
the critical exponents near the demixing point of directed polymer blends and compare 
them to those for point particles. Analogous phenomena should occur in binary mixtures 
of two different superfluids. Among the possible candidates would be films of ffe 6 and 
He 4 or possibly, a superfluid monolayer of spin aligned hydrogen mixed with superfluid 
He 4 . As the mass ratio, mi/m2, of the two species is increased towards unity with h > 0, 
quantum fluctuations should eventually cause the two superfluids to mix. The same effect 
may be achievable by varying the properties of the substrate or formally, by increasing %. 
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A phase diagram as a function of temperature, composition and ft is shown in Figure 4. 
Note the incomplete phase separation appearing for nonzero ft. The theory applies to the 
critical point at T = 0. In Section 5 we will analyze such binary mixtures in greater detail. 
Although we shall concentrate on 2 + 1 quantum systems, our results are easily adapted 
to 3 + 1 dimensional quantum problems which are well described by mean field theory (see 
Section 4). 

1.5. Outline 

In Section 2 we formulate the Landau theory that will be the focus of this paper. 
We follow the derivation developed in for flux lines, which leads to a coherent state 
functional integral. From there we describe the hydrodynamic theory with and without 
free ends as in H. In appendix A we present a similar derivation for two interacting 
polymer species. 

In Section 3 we describe the correlation functions of this theory which could be mea- 
sured by X-ray or neutron diffraction experiments. We also discuss the meaning of the 
direction dependent compressibility which characterizes these directed line systems. 

In Section 4 we study the critical behavior of our model via an expansion in 2 — e + 1 
dimensions. We reproduce the results of Larkin and Khmel'nitskii for e = when there are 
no free ends @, and find that free ends are a relevant perturbation leading to the critical 
behavior of a (d + l)-dimensional Ising model with short range interactions. 

Finally, in Section 5 we discuss the critical mixing of two polymer or quantum boson 
species. Using the free energy derived in appendix A, we show that, just as in the point- 
particle case, the critical mixing of line-like objects is in the same universality class as 
single component critical phenomena. 



2. Coherent State Functional Integrals and Landau Theory 

The derivation of (|1 . 1] ) from a more microscopic free energy for directed lines follows 
the presentations in H and |1. We only outline the main points here, first reviewing the 
theory for chains which span the system. Internal free ends are then introduced by adding 
a source to the boson coherent state field theory 0,[| . 



We start from the following path integral partition function for the N chains, 



'N 






exp 



"-B-* JQ 



dz 



N 

E 



:§(£)' + |>'M -*«']} 



(2.1) 



where we have restricted our attention to interactions between lines through an equal 
"time" potential V[rj(z) — Tj(z)]. 

The boson field theory is derived in three steps. One first writes the transfer matrix 
associated with ( |2.1[ ) in terms of a many-particle real-time Schrodinger equation for the 
polymer positions, and then second quantizes this equation using bosons. Finally, we pass 
to a coherent state functional integral representation of the second quantized formalism. 

The grand canonical partition sum which results is 



oo „ 

Z gr = J2 e L ^ N/kBT Z N = / [#][#* 

7V=0 J 



-S[VM 



(2.2) 



where i/j(r, z) is a complex boson field. The boson action S reads 



S = I dz I d z r 



r(r,z)(d z -DV ± 2 -fi)ij(r,z) + 
- / d 2 r' v(r — r')\ip(r, z)\ 2 \ip(r', z)\" 



(2.3) 



where D = k s T/2g, \x = n/k B T and v = V/k B T. Standard manipulations show that the 
density of flux lines is 

p(r,z) = mr,z)\ 2 . (2.4) 



2. 1 . Entangled Limit with Free Ends 

We now assume that all phases are dense enough to be entangled and expand around 

l —1/2 

the mean density po • For polymers of length £, entanglement means yDl >> p , 
while for quantum systems, we require very low temperatures so that the thermal de 
Broglie wavelength is large compared to the particle spacing. We parameterize the field if) 
by i> = \fp^ i a change of variables with unit functional determinant, and find from 
that, 



S 



dzd r 



d(Vlp) 



+ D P (yj) 2 + i P d z e -pp + -p 2 

4p 2 



(2.5) 



Upon introducing an auxiliary field t(r,z) conjugate to Vj_0(r, z) via a Hubbard- 
Stratonovich transformation, and integrating out 6(r,z), we arrive at the partition sum 
( |1.2|) P,|2l[|. See Appendix A for the analogous transformation applied to binary mixtures 



of quantum bosons. 

We can add a (Poisson distribution) of free ends to the theory by adding to the action 
( P75| ) a term of the form || 

(2.6) 
dzd 2 r Irj^fp cos 9 



The mean polymer length £ is related to r\ by 2rj£k B T = y/po [Q. We may expand this 
term in powers of 9 if the fluctuations in 9 are small, which they will be if the polymers 
are sufficiently dense and entangled. To leading order in 9, equation (|2.6|) becomes 



8S « const. + / dzd 2 rr]^/p9 2 (2.7) 

The Hubbard- Stratonovich transformation now leads to 

S = /^{f + ^-^ + ^ + V,t]*} (2.8) 

Upon expanding around an average density po, we arrive at the theory for finite length 
polymers (|1.8|), with G = \j4r\Jp~Q and h = 2D/p . 



3. Correlations and Compressibility 

From the free energy ( |1.8| ) we can, by considering only quadratic terms, easily find the 
hydrodynamic form of the structure function which could be measured by either neutron 
or x-ray diffraction. Here and henceforth we change units so that e = 1. The Fourier- 
transformed density-density correlation function is then 

S{q±,q z ) = 2 

P ° (3 1) 

KT{h + Gql) K ' J 



(Gaql + Gh + ah) q 2 z + (Gb + h) q\ + Gqi + hb ' 

Note that G drops out of ( |3.1|) upon setting h = 0. This would have to be the case in 
the boson analogy, since turning h = neglects quantum mechanics by decoupling the 



momentum from the position. The constraint (|1.3| ) is unimportant in this classical limit. 



In the limit G — ► oo , we recover the correlations typical of infinite polymers with no 
heads or tails from (|3~T|), 

S{Q± ' ^ = (h + aqiui+bql+qf ^ 

For small momenta, we can examine the curves of constant S(q ± ,q z ). It is easy to see 
that these curves will go into the origin on straight lines, i.e. q z ~ q ± . This is a generic 
feature expected to hold for all systems with the morphology of directed polymers, where 
the direction is explicitly chosen by an external field. Recalling that the limit as \q\ —>■ of 
the structure function is related to the compressibility of the sample |2!| , we can determine 



the compressibility of a sample of directed lines, which will depend on the direction of the 
compression ||23|1 . Upon considering a compressional wave at an angle <fi from the z axis, 



we take the momentum to be q z = q cos (p and q± = q sin <fi, and find for small q, 

k T 

S(q±,q z ) ~ -7 -2 = k B T X T (0), (3.3) 

i^° b + h (cot (/)) 

where Xt is the compressibility, and is a function of the direction <p in which the compression 

is made. Note again that the h = limit is especially simple: The compressibility no longer 

depends on the compression angle if tangent fluctuations are ignored. We then recover the 

result for an isotropic fluid of point particles, namely Xt = b~ x . The compressibility enters 

the linear response relation: 

(Sp(q±, Qz)} = X(q±, q z )Sa(q ± , q z ), (3.4) 

where da is a sinusoidally varying force directed along q = (q±,q z ). The unusual (jy 



dependence of this response function may be understood as follows: Equations (^J) and 
( P73| ) predict that (dp) = for compressions when <fi — > 0. This result arises because 
squeezing the system along the z-axis will not change the number of polymers piercing 
any constant z-slice. In general, the linear response ( |3.4|) depends on the directions of the 
compression relative to z, as parameterized by 4>. 

Returning to finite length polymers, and hence to finite G, we repeat the above analysis 
and find a ^-independent result 

k T 
S(q±,q z ) ~ ^— • (3.5) 

Because at sufficiently large distances (larger than the average polymer length £) finite 
length polymers will ultimately behave like point particles, the long-wavelength properties 
of the sample will not exhibit any character of the directed constituents. At intermediate 
scales, the factors of 0{Gq\) may not be neglected in (|3.1|), and we return to the direction 
dependent result (|3.3|). 
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4. Renormalization Group and Finite Length Polymers 

In the last section we saw that at long- wavelengths the line-like nature of the directed 
melt disappears for finite chain lengths. Here we construct a renormalization group which 
allows us to study this crossover in detail. When G — » oo, we recover the results of Larkin 
and Khmel'nitskh 0. For G finite, however, the system crosses over to the behavior of 
a (d + 1) -dimensional Ising model with short range interactions [[20]. We shall treat the 
^-direction separately from the directions in the plane perpendicular to z. We let d be 
the dimension of the bosons so that the corresponding directed polymers exist in (d + 1)- 
dimensions (see Figure 1). 

Consider first the graph which renormalizes the four-point interaction dp 4 . To inter- 
pret this graph we need to invert a (d + 1) x (d + 1) matrix to get the propagator, as 
there are d degrees of freedom in t and an additional degree of freedom from the density 
fluctuations, dp. However, as mentioned earlier, the transverse modes of t decouple. To 
make this explicit, we decompose t in momentum space, 



t{q±,q z ) = t L + t T , (4.1 

with 

tL.i = - 

q. 

and 



t L j = ^rt k (4.2) 



trj = [S jk - ^J t k (4.3) 

In (4.2) and ( f4.3|) the indices j and k only run over the <i-dimensional space perpendicular to 



z and we use the summation convention. Additionally we take t^j = —iqj^(q±, q z )- With 
this parameterization, the fourier transform of the divergence of t is simply iq ± ■ t — q\n, 

and the fourier transform of the t 2 term in ( |1.8|) is J dq z d d q ± [q^n 2 + tx 2 h-(q± • tr) 2 ]- 

Thus we can integrate out tj- from ( |1.8| ), and we are left with a theory only involving dp 
and 7r. 

Inverting the 2x2 quadratic form which remains gives the following useful propaga- 
tors: 

(\5p(q^q z )\ 2 ) = ( Gag 2 +Gh + ah) q 2 + (G J + h )ql + Gqi + hb ^ Aa) 

(U( Q a)\ 2 ) = - gj + 6 + aql + Gq\ 

M Wi, qz)\ 1 g 2 (Gag 2 +Gh + ah ) q 2 + ( G6 + h)ql+ Gq 4 +hb { ] 

(6 P *(q ± ,q z Mq^q z )) = ^^ + - - — } ^ + ^ - - } g + ^ + - (4.4c) 
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After inserting these propagators into the graph shown in Figure 5, we integrate out the 
frequency- like variable q z and find the renormalized four point interaction associated with 
equation ( |1.8|) is c R = c + 5c + C(c 3 ), where 



5c 



3c 2 



8 r(f)(47r) d / 2 



- X dq, 



1 + 



h 



Gq< 




(4.5) 



Near the critical point, b m and so when G = oo, 5c diverges in the infrared for d < 2. 
This infrared divergence suggests an e-expansion around d — 2, where e — 2 — d. Our 



approach follows |24[] , which is a variation of the dynamic renormalization group of [55] . 
First we integrate out a momentum shell in q ± from Ae~ £ to A, but integrate freely over q z . 
We then rescale our variables so that the ultraviolet cutoff in the perpendicular direction 
is held fixed. After rescaling p and t accordingly, we are left with the same theory but 
with different coupling constants. When this procedure is iterated, c is driven toward a 
fixed point which describes the universal long wavelength critical behavior for large G. 

4-1. Momentum Shell Integration 

We must first integrate out the transverse momentum in the range Ae~ £ < q ± < A. 
This can be done straightforwardly by expanding the functional integral in c. Care must 
be taken to account for all possible contractions of the operators in the expansion. The 
symmetry factors can be found in the usual way for Wick expansions. It is important to 
note that diagrams renormalize the remaining low momentum modes, and are not simply 
expectation values. Graphs which renormalize the quadratic and quartic contributions to 
Landau theory are shown in Figure 6. Upon carrying out this procedure, we arrive at the 
following relations for the intermediate values of the coupling constants in ( |1.8| ) 

ti = h 

a' = a 



b' = b + 






V 






(4.6a) 
(4.66) 


c r A d d q ± 
4VG Ja^ (27r) d 


h + Gq 2 x 1 


(4.6c) 


V 


/b + ql ^ h + aq 2 ±+a h^ 


\ 


, 3c r A d d q ± 
SVG J A e-' (^) d 


yfh + Gq 2 ± 1 




] (4.6d) 


V( b + <ll) 3 ^h + aq 2 ± + 


«B. 


(4.6e) 
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where we have gone to one-loop order. In Q4.6|) we have already integrated over q z . 

We rescale perpendicular lengths by R' ± = R ± e~ e and the z direction by Z' = 
Ze Jo 1( ' . In momentum space, these rescalings read q' ± = q ± e^ and q' z = q z e<>o 7l , 
where the function j(£) is to be determined. 

When doing the first momentum shell integration, the coupling constants were inde- 
pendent of length scale. However, they then acquire a momentum dependence because we 
have absorbed the large momentum effects into them. The correct renormalized theory 
is a coupled set of integral equations where the coupling constants are taken to be scale 
dependent. An alternative, but equivalent approach is to integrate over a small momen- 
tum range where the coupling constants are approximately fixed and then repeat the entire 
calculation iteratively. This leads to the usual differential renormalization group equations. 

4-2. Recursion Relations and Critical Behavior 

We now choose an infinitesimal momentum shell e _(5 , and take the limit 5^0. This 
leads to differential renormalization group equations which can be integrated to produce 
the couplings appropriate for a cutoff Ae~ e in the perpendicular direction. We use units 
such that A = 1 in the following. We also define the renormalized "reduced temperature" 

r = b+ 2vh; A - (4 ' 7) 

which vanishes at the critical point in the fluctuation-corrected theory. The constant 
Ad = 2/[r(|)(47r) d ' 2 ] is a geometrical factor. The differential recursion relations valid 
near d = 2 are 



dh( - £) --h(4-2y) (4.8a 



d£ 
da(£) 



a (2- 2 7 ) (4.86 



drU) l c \ ch 

_JLZ = r 2 + (4.8c 

ai \ 16n^/h + a + a^J 16nG^/h + a + a^ 

M \ IQiiJh + a + a^ 

^f = G(2-2 7 ) ' (4. 8e 
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We set 7 = 2 in order to hold h fixed, and note that a is an irrelevant variable. 
We have expanded about G _1 = 0. Examining these recursion relations, we see that 
perturbation theory amounts to an expansion in two parameters, namely c = c/\/h + ai± 
and (G) -1 = hG~ l . These couplings have the recursion relations 

where e — 2 — d. 

We first consider the subspace of theories with G = oo. In this case there is a stable, 
nontrivial fixed point for positive e. Though we have already noted that this theory is 
identical to the uniaxial ferroelectric, and hence should have all the same critical behavior, 
a direct calculation of these exponents provides a useful check of our renormalization group 
method. In particular we can calculate the logarithmic corrections to the compressibility 
discussed in Section 3 in the critical dimension of 2 + 1 = 3. Upon setting e = in ( |4.9q|) 
we solve for c(£) and find that 

m = ~%r; <4 ' 10) 



which leads via (|4.8dj to 



Mi, , 35 (°) 



-1/3 

rW = r(0)e 2 Ml + ^l (4.11) 

Since at this order in e there is no nontrivial rescaling of the fields, we have 
Xfex, q z ; ro) = / dzd d x e^'*-*'* (p(x, z)p(0, 0)) e=0 

= e" Xl (e^Jo^ M ' qz]r (i) 



(4.12) 



We now choose £ = £* large enough so that r(£*) has grown to 0(1). Then we have from 

1 = r {t) = r (3c /167r)- 1/3 e 2r {£* )" 1/3 (4.13) 
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so that 



and thus 



^ ^-MwjVs 



~ r. 



o _1 |lnr | 1/3 (4.14) 



X(q sin <p, q cos </>; r ) ~^ : r7r; - ; — — —% (4.16) 



I lnrol 1 / 3 
X(q±,(lz; ro) = X(^x9x,6?«; !) ( 4 - 15 ) 

where £ 2 ~ £ 2 ~ (lnro) 1 / 3 /^, in agreement with the results of 0. Note the anisotropi- 
cally diverging correlation lengths produced by the directed line-like nature of the degrees 
of freedom. The two-loop calculation of [^6[ leads to the result £ z ~ £ x 243 to 0(e 2 ). 
Recalling the discussion in the last section of the compressibility, we have finally 

lnrol 1 / 3 1 

9^0 r l + Zidlnrol 1 / 3 /^) (cot. 

which should be compared to the mean field result (|3.3|). Note that the effective direction 
of the compressional wave changes as we go to longer length scales. Near the critical point 
we always find x(</>) ~ h( TT& i which appears to be a universal result. 

Similarly, in the ordered phase, we can calculate the specific heat exponent by flowing 
along a renormalization group trajectory and matching on to an effective free energy at 
long-wavelengths (A more complicated calculation shows that this behavior occurs above 
T c as well). In the ordered phase, the bulk free energy is approximately — 3Q(r 2 /2c), where 
O is the total volume of the system. We then have 

r(£) 2 



(4.17) 



Fe = — 30/ 

2c(f) 

2 c(£) 

Again choosing £ = £* as in ( |4.13| ) we find that 

Ft = -\e- M \~ C 7( ^'0o^y ~ e-V+Wt ~ r 2 | lnr \^ (4.18) 

Differentiating twice with respect to ro, we recover the correct specific heat-like behavior 
at the critical point. 

Going back to the full theory we now allow G~ l to take on finite values. In this 
case we are driven away from our fixed point and G~ x flows towards oo. This reflects 
a crossover to an anisotropic but point-like phase where the constraint ( p..3|) , and hence 
the line-like nature of the polymers is unimportant. Because t decouples from the density 
fluctuations the quantum mechanical bosons behave classically at long wavelengths. A 
similar crossover occurs in polymer nematics when hairpins are introduced into the theory 
and the previously directed polymers become isotropic J|. 
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5. Mixing Exponents 

We now consider directed polymer blends, or, equivalently, binary mixtures of two 
superfluids. As with binary mixing in classical systems of point particles, the mixing 
exponents are identical to the liquid-gas exponents. We demonstrate this equivalence by 
showing that the Landau- Ginzburg theory controlling the mixing fraction is the same as 
that for a single component critical point, (|1.1|) . For simplicity, we restrict our attention 
to the limit G — > oo. In appendix A we derive the quadratic part of the free energy of 
two polymer species using the boson representation. To these terms we add nonlinear 
couplings to account for the configurational entropy of each polymer species, multiparticle 
interactions, etc., and consider the model free energy 



dzd d x 



hi f 2 
2 



^tl + hl 



*1 + "^2 + fcl2tl-t 2 + ^(VJpi, 



1 

2 l 



+ o(Vx5p 2 ) 2 



+ j Spj + —6pl + b 12 5 Pl 5p 2 + —Spl + -^5p% 

. C\ 4 C 2 - 4 Cl2 2r 2 

+ ^j-opi + 4j-5p 2 + — S Pl 5p 2 



(5.1) 



subject to the constraints 



djpj + Vx-tj 







J = 1,2. 



(5.2) 



We have neglected gradient terms of the form (d z 5pi) 2 , because these turn out to be 
irrelevant at the demixing point. Terms proportional to 5p±5p 2 and 5p 2 5p\ could have 
been added, but do not affector our results in any essential way. Additionally, we assume 
that the (V x 5pi) 2 terms have been diagonalized so that the term Vx^pi • V ± 5p 2 is not 
present. It is useful to pass to sum and difference variables, i.e. 5p ± = (5p\ ± 5p 2 )/\/2. 
We only keep terms in the free energy up to second order in the noncritical mode 8p + , but 
keep all terms up to quartic order in 5p_. We also define fields t ± = (ti ± t 2 )/v / 2 which 
are the tangents associated with the sum and difference variables. The free energy now 
takes the form 



dzd d x 



h, 9 

— tl + 
2 + 

b + 



h. 



t 2 + h 



+ _t + t. 



+ ^(Vx5p H 



+ ^(Vx5/9_ 



!■ f 5 pl + Y sp2 - + b+ ~ 6p+5p - + lf 5/0 - + % 6p4 



(5.3) 



with the constraints now reading 



d z 5p ± + V ± -t ± = 0. 



(5.4) 
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We can again integrate out the transverse parts of t ± , and use the constraint (|5.4| ) to 
solve for the longitudinal parts. The resulting quadratic part of the free energy reads (in 

momentum space) 

n / J \ 



quad 



dq z d d q ± 

2n (2n) d 



U h 4 +qi + b < 



1 



<ll 



+ ^[h.^+qi+b\\5p 



2 V q 

+ h. 



£ 



5p + 5p„ + b + _5p + 5p_ 



(5.5) 



At this point we are in the position to integrate out Sp + , resulting in 



dq z d d q ± 

2tt (2nY 



h_^ + q 2 ± +b_ 






h^ + ql+b, 



+ 2 



\8p. 



(5.6) 



Upon expanding the denominator of ( |5.6|) in q ± and q z , keeping only relevant and 
marginally relevant terms, and reintroducing t_ we finally have 

1 (. 2b+_h+_ b 2 , h, 

1- + ~ I t~ + 

2 



"J 

subject to 



dzd d x 



h 



K 



t 2 + - i + 



b\ 



(VJp- 



1 U b l 

+ 2 [ b -'t 



6£ + ^ + c ^ P t 



(5.7) 



(5.8) 



d z 5p_ + V ± t_ =0. 

As in the single component model, the critical mixing line is gotten by setting w_ = 0. 
This could be done by altering the chemical potentials of the two species, for instance. 
Our model is thus equivalent to the original model ( p. . 1|) , except for the unimportant 
coupling a(d z Sp_) 2 . We conclude that the critical mixing of line-like objects is also in the 
universality class of the three-dimensional uniaxial ferroelectric. 
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Appendix A. Derivation of Hydrodynamics for Binary Mixtures 



Our starting point is a pair of boson field theories, each describing a distinct superfiuid 
or polymer species @. We then couple them through their densities, as well as their 
"currents" or tangent fields. More precisely the partition function is 



[#!*][#!] [# 2 *][# 2 ]e- 



-S[lp*,1pl,1pZ,1p2 



where the action is 
F 



The individual terms read 



s^s^^H + s^lH + SubPl^i^lH- 



Sj = I dzd d x 



B, 



^(d,-z^v x 2 -^)^ + ^h/V 4 



while the interaction is 

dzd d x [S^I^H^I 2 + ^712 (</>i* V x ^i - V'lV^i*) • (^v^ 2 - ^v^)] 



The change of variables 



fa 



'Pi* 



leads to 



S 



J dzd * X < 



4pi 4 

ip2d z 9 2 + D2{y±p2) + D 2 p 2 {VJ 2 ) 2 - U2P2 + ^p\ 
4p 2 4 

k + S12P1P2 + r] 12 pip 2 V ± 9i ■ V ± 6 2 



(A.1) 



(A.2) 



(A.3) 



(A.4) 



(A.5) 



(A.6) 



We now introduce two fields which we shall see are proportional to the tangent fields of 
the individual species. We start with the identity 



expUdzd d x [D lPl (VJ 1 ) 2 +D 2 p 2 {VJ 2 ) 2 + T ll2 p 1 p 2 VJ 1 ■ V ± 9 2 ] 

{[" D2P2PI + -D1P1P2 ~ 2?7i2PiP2PrP2 
J dzd d x D 1 D 2PlP2 - rfi 2 p\pl 

I + ip 1 -V ± 8 1 + ip 2 -V ± 8 2 



(A.7) 
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Upon using this identity to replace the terms quadratic in 6 appearing in ( |A.6|) and inte- 
grating over 6j, we find that 

Z = J[d Pl }[dp 2 ][dp 1 }[dp 2 }5[d zPl + V ± - Pl ]5[d zP2 + V x - P2 ]e- s '^^' Pl ' P2 ] (A.8) 



where 



( D 2 p2V 2 i + ^lPiPi - 277 12 piP2pi-p2 



5' 



2 ^2^2 



dzd x < 



D l D 2 pip 2 - VI2P1P2 

| Dl(VxPl) 2 | ^ 2 (VxP2) 



4/0i 



4p 2 



MlPl - ^2P2 + ~^p\ + ~^p\ + B 12 p!p 2 



(A.9) 



We now expand pi = pcM + dpi around the minimum of the bulk free energy and rescale 
dpj so that the kinetic term in the action is 2} } T (V ± Spj) 2 to obtain 



S' 



KT 



dzd d x I 2 



■7ft? + -#*! + ^i2ti-t 2 + J(V x 5pi) 2 + ^-(V ± 5p 2 ) 2 



1 
2 



1 

2 l 



Ol 9 °2 c 2 1 r r 

+ TT°Pl + TT"/°2 + V12OP1OP2 



(A.10) 



where t j = ^Dj/2pj^ pj and the /ij and 6j can be read off from ( |A.9| ) in the expansion of 
Pj around pj^. The delta-functionals in ( |A.8| ) impose the constraints 



d z Spj + V x -tj = 0. 



(A.ll) 
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Figure Captions 

Fig. 1. Hydro dynamic volume averaging over directed lines in (d + 1) dimensions which 
leads to the coarse-grained density and tangent fields used in this paper. 

Fig. 2. Speculative phase diagram for Type II superconductors with an attractive inter- 
action between flux lines which allows for coexisting dense and dilute entangled 
flux liquids. At sufficiently high temperatures, this line terminates in a critical 
point. The critical isochore, which marks the extension of this line into the single 
phase region, is shown as a dotted line. 

Fig. 3. Pressure-temperature phase diagram for a quantum fluid with the unusual pair 
potential shown in the inset. The outer minimum dominates in the superfluid 
phase Si, while the inner minimum is primarily occupied in the higher pressure 
superfluid S2. Phases 5*1 and S2 are separated from a normal liquid N by the 
dashed A-line. Points a and b are conventional critical points. By enhancing 
quantum fluctuations at T = (say, by increasing fi), on arrives at the critical 
point c discussed in the paper. 

Fig. 4. Phase diagram for a binary mixture of two superfluid species as a function of 
the concentration of species B, temperature and Planck's constant. This paper 
addresses the nature of the critical point which occurs at T = 0. 

Fig. 5. Hartree graph which diverges for d < 2. We integrate here over all internal 
momenta. 

Fig. 6. Graphs which contribute to the renormalization group calculation. In these 
graphs, the internal lines are integrated over all values of q z but q ± is only inte- 
grated in a momentum shell between Ae~ e and A. Figure (a) is a graph which 
renormalizes the propagator, and figure (b) is a graph which renormalizes the 
four-point coupling c. 
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